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^— s ^ We investigate the infrared behavior of the spectrum of scalar-dressed, asymptotically Anti de 

fvj . Sitter (AdS) black brane (BB) solutions of effective holographic models. These solutions describe 

scalar condensates in the dual field theories. We show that for zero charge density the ground state 

of these BBs must be degenerate with the AdS vacuum, must satisfy conformal boundary conditions 

for the scalar field and it is isolated from the continuous part of the spectrum. When a finite charge 

density is switched on, the ground state is not anymore isolated and the degeneracy is removed. 

Depending on the coupling functions, the new ground state may possibly be energetically preferred 

with respect to the extremal Reissner-Nordstrom AdS BB. We derive several properties of BBs near 

extremality and at finite temperature. As a check and illustration of our results we derive and 

discuss several analytic and numerical, BB solutions of Einstein-scalar-Ma^cwell AdS gravity with 

different coupling functions and different potentials. We also discuss how our results can be used for 

understanding holographic quantum critical points, in particular their stability and the associated 

t*H' quantum phase transitions leading to superconductivity or hyperscaling violation. 
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I. INTRODUCTION 

Holographic models have been widely used as a powerful tool to describe the strongly coupled 
regime of quantum field theory (QFT) |ll428|. In the usual holographic setup, the 4-dimensional 
gravitational bulk is described by an asymptotically anti-de Sitter (AdS) black brane (BB), gener- 
ically sourced by a (neutral or charged) scalar field and by an electromagnetic (EM) field. The 
usual rules of the AdS/CFT correspondence are then used to describe a dual 3-dimensional QFT 
at finite charge (or finite chemical potential) and where some scalar operator may acquire a non- 
vanishing expectation value. In particular, when the BB solution is sourced by a nontrivial scalar 
field configuration in the bulk, the dual QFT shows a (neutral or charged) scalar condensation. 

These effective holographic theories (EHTs) have a wide range of applications. They have been 
used to give a holographic description of many interesting quantum phase transitions, such as 
those leading to critical superconductivity and hyperscaling violation [l|, 0, mM, Ej Hi ^M, |30| ■ 
The properties of the QFT at small temperatures are essentially determined by the quantum 
phase transition at zero temperature. Thus, understanding the quantum critical point provides a 
full characterization of the thermodynamical phase transition at small temperatures. 

If one assumes, as we do in this paper, that the asymptotic geometry is the AdS spacetime, the 
dual QFT shows a universal conformal fixed point in the ultraviolet (UV). The nontrivial dynamics 
therefore occurs in the infrared (IR) region, at the corresponding critical points. In a Wilsonian 
approach, EHTs should be first classified in terms of fiows, driven by relevant operators, between 
critical points corresponding to scale-invariant (more generally scale-covariant) QFTs. Two other 
relevant characterizations of the critical points are: a) the distinction between fractionalized phases 
(sourced by non-zero electric flux in the IR) and cohesive phases (sourced by zero electric flux in 
the IR); b) phases with broken and unbroken C/(l) symmetry [la. [28|. 

Progress in the classification and understanding of IR critical points have been achieved following 
various directions. In particular, it has been shown that in the case of hyperscaling preserving and 
hyperscaling violating solutions, quantum critical theories may appear as fixed lines rather than 
fixed points j28| . Hyperscaling preserving solutions appear indeed as fixed points and correspond 
to AdS4, AdS2 X R^ and Lifshitz bulk geometries. However, hyperscaling violating solutions are 
characterized by an explicit scale and therefore appear rather as critical lines generated by changing 
that scale or equivalently the charge density [1^ [28l - t30fl . 

A crucial point for understanding these quantum critical points is the presence of a scalar 
condensate. Indeed nontrivial configurations of (generically charged or uncharged) scalar fields 
play several crucial roles: (i) nontrivial scalar fields are dual to relevant operators that drive 
the renormalization group (RG) flow from the UV flxed point to the IR critical point (or line); 



(ii) scalar fields are the sources that sup por t the IR hyperscaling violating geometry allowing for 
both fractionalized and cohesive phases jlj!, [2^; (iii) charged scalar condensates break the C/(l) 
symmetry and generate a superconducting phase |16l |28| . 

Despite the recognized relevance of scalar condensates for describing holographic critical points, 
we are far from having a complete understanding of the physics behind them, in particular we have 
very few information about their stability. For instance, one would like to understand why at zero 
(and small) temperature the hyperscaling violating phase is energetically preferred with respect 
to the hyperscaling preserving phase. In this paper we will move a step forward in this direction 
by asking ourselves a simple, but relevant question: what is the energy of the ground state of a 
neutral asymptotically AdS BB sourced by a generic nontrivial scalar field? We show that for zero 
charge density the BB ground state must be degenerate with the AdS vacuum. This degeneracy is 
the result of an exact cancellation between a positive gravitational contribution to the energy and 
a negative contribution due to the scalar condensate. 

Moreover, we also show that for the BB ground state the symmetries of the field equations 
force conformal boundary conditions for the scalar field, i.e. boundary conditions preserving the 
asymptotic symmetry group of the AdS spacetime. The conformal boundary conditions correspond 
to dual multitrace scalar operators driving the dynamics from the UV conformal fixed point to 
the IR critical point. In the case of an IR hyperscaling violating geometry sourced by a pure 
scalar field with a potential behaving exponentially, a scale is generated in the IR. On the other 
hand we will show that, in the case of pure Einstein-scalar gravity at finite temperature T , the 
boundary conditions for the scalar are determined by the dynamics and are, therefore, generically 
nonconformal. This means that the ground state for scalar BBs is isolated, i.e. it cannot be 
obtained as the T ^ limit of finite-T BBs with conformal boundary conditions for the scalar 
field. 

When a finite charge density p is switched on, the degeneracy of the ground state is removed. 
Because an additional degree of freedom (the EM potential) is present, the boundary conditions 
for the scalar field are not anymore determined by the dynamics. The freedom in choosing the 
boundary conditions arbitrarily can be used to impose conformal boundary conditions also for BBs 
at finite temperature. The ground state for scalar BBs is therefore not anymore isolated from the 
continuous part of the spectrum. The coupling between the bulk scalar and EM field determines 
if it is energetically preferred with respect to the extremal Reissner- Nordstrom (RN) AdS BB. 

We also derive several properties of the scalar BB near extremality and at finite temperature. 
For instance, we show that scalar-dressed, neutral (charged), BB solutions of radius r^ (and charge 
density p) only exist for a temperature T bigger than the temperature of the Schwarzschild-AdS 
(Reissner-Nordstrom AdS) BB with the same r^ (and with the same p). 

As a check and illustration of our results we give and discuss - both analytically and numerically - 
several (un)charged, scalar-dressed BB solutions of Einstein-scalar-Maxwell AdS (ESM-AdS) grav- 
ity with minimal, nonminimal and covariant coupling functions and different potentials (quadratic, 
quartic, exponential). 

Finally, we also discuss the relevance of our results for understanding holographic quantum 
critical points, in particular their stability and the associated quantum phase transitions. 

The structure of this paper is the following. In Sect. [IT] we present the general form of the EHTs 
we consider. In Sect. IIIII we investigate the spectrum of this class of theories in the IR region. In 
Sect. IIVI we derive extremal, near-extremal and finite-temperature BB solutions of pure Einstein- 
scalar gravity theories in the case of a quadratic, quartic and exponential potential. We also derive 
their thermodynamical behavior and their critical exponents. In Sect. |V] we derive and discuss 
charged solutions with the scalar minimally, nonminimally and covariantly coupled to the EM field. 
Finally in Sect. IVII we end the paper with some concluding remarks about the relevance that our 
results have for understanding the dual QFT, holographic quantum critical points, in particular the 
stability of the latter and the associated quantum phase transitions leading to superconductivity 
or hyperscaling violation. In Appendix \X\ we discuss perturbative solutions in the small scalar field 
limit. 



II. EFFECTIVE HOLOGRAPHIC THEORIES 

We consider Einstein gravity coupled to a real scalar field and to an EM field in four dimensions: 



d'^x 



R-l{d<f>f 



z{4>) 



F^ - ¥{(/)) - Y{(I))A'^ 



(2.1) 



where Ff^^^ = d^A^ — dyA^ is the Maxwell field-strength. The model is parametrized by the 
gauge coupling function 2'(0), by the self-interaction potential V{(t>) for the scalar field and by the 
coupling function Y{(f)) giving mass to the Maxwell field. 

The action (|2.ip defines ESM theories of gravity, which are also called EHTs and are relevant 
for holographic applications. EHTs have been widely used to give a holographical description of 
strongly-coupled QFTs with rich phenomenolo gy s uch as quantum phase transitions, superconduc- 
tivity and hyperscaling violation |ll. I2I. la-lol. Ha. IiqI. l28l - (30| . Moreover, models like (|2.ip generically 
appear, after dimensional reduction, as low-energy effective string theories. The action (|2.ip can 
be also interpreted as an EHT for a complex scalar field that enjoys a C/(l) symmetry [23 ■ In this 
context the real scalar describes the modulus of the charged scalar and the phase with broken 
(unbroken) f/(l) symmetry is obtained by F ^ (y = 0). 

Although our considerations can be easily extended to the case F 7^ 0, we will focus for simplicity 
on the case of unbroken U{1) symmetry, Y — Q. We will briefly comment on the case y 7^ in 
Section Eel 

We are interested in electrically charged BB solutions of the theory, i.e. static solutions with 
radial symmetry for which the topology of the transverse space is planar. Using the following 
parametrization for the metric: 



ds' 



dr'^ 
-X{r)dt'+ +H\r){dx'+dy'), 

A(r) 



(2.2) 



the Einstein and scalar equations read: 



El = JA 

H 

{XHH'Y = -H 



-, {XH 



2\M 



V ZA'^'' 



<-v,'^HM;f-ff 



-2H^V, 



(2.3) 
(2.4) 

(2.5) 



The ansatz 
solved for A'^. 



is very convenient, as in these coordinates Maxwell's equations can be directly 



A' - ^ 



(2.6) 



where p is the charge density of the solution. Note that only Eqs. (|2.4[) and (|2.5p depend on the 
EM field and only through A'q. Therefore, substituting the solution above into the remaining field 
equations, we can completely eliminate the EM field and solve Eqs. (|2.3p - (|2.5p for A, H and (j>. 

We will consider models for which the potential F (0) has a maximum at = and Z'{(p = 0) = 0, 
with the local mass of the scalar ml = V"{0) satisfying the condition rri^p < m^ < — 2/L^ and 
with V{0) = — 6/i^, where rrigp ~ — 9/(4i^) is the Breitenlohner-Freedman (BE) bound [3l| in 
four dimensions and L is the AdS length^. The presence of an extremum of V{(l)) and Z{(l)) at 



^ The results of our paper can be easily extended to the scalar-mass range m^p < 
dual CFT to is known to be unitary. 



mi < mi 



1/L^, where the 



= implies the existence of a Reissiier-Nordstrom-AdS (RN-AdS) BB solution, 



M p^ 
12 2^ + 4^2 



^=72-^ + 73' H^r, = 0, (2.7) 



which is characterized by a trivial scalar field configuration. 

On the other hand well-known "no-hair" theorems [32h35| make the existence of BB solutions 
endowed with a nontrivial scalar field a rather involved question. Scalar-dressed solutions are 
particularly important for holographic applications because they describe dual QFTs with a scalar 
condensate. 

The AdS, r = cx), asymptotic behavior requires the following leading behavior of the metric and 
the scalar field: 

ds^ = --p^df + -^dr"^ + r^{dx^ + dy'^) 



L2 



r 






with Ai_2 = — — — 5 — ^ — . Because the AdS spacetime is not globally hyperbolic, this asymptotic 
behavior must be supported by boundary conditions on Oi, 02- Dirichlet boundary conditions 
that preserve the asymptotic isometrics of the AdS spacetime are Oi = 0. However, in the range 
of scalar masses rn%p < ml < m?gp + l/i^, boundary conditions of the form (/ is a constant): 

Ol = /o^/'^^ (2.9) 

which preserve the conformal, asymptotic symmetries of the AdS background, are also allowed |36l |. 
More in general, boundary conditions of the form 

Oi^W{02), (2.10) 

can be used. For a generic form of the function W the asymptotic AdS isometrics are broken, yet 
an asymptotic time- like Killing vector exists and the gravitational theory admits a dual description 
in terms of multitrace deformations of CFTs |37H43| . 

Apart from their UV AdS behavior, the scalar-dressed solutions of EHTs are also characterized by 
their, small r, IR behavior. This IR behavior is of crucial relevance for holographic applications, in 
particular in the context of the AdS/Condensed Matter correspondence [5|. Generically, we expect 
the IR regime not to be universal, but rather determined by the infrared behavior of the potential 
V{(t)) and of the gauge coupling functions Z(0), Y{(f)). Nevertheless, we will discover in the next 
sections some features of the IR spectrum of EHTs, which are model-independent and related to 
the scaling symmetries of the UV AdS vacuum. 

Although we will be concerned with general features of EHTs, for the sake of definiteness we will 
mainly focus on three classes of models with different IR behavior of the potential V{(f)): 



a) The potential has a quadratic form 






V{c^) = -i2+^- (2-11) 



This corresponds to the simplest choice for the potential, which has been widely used in 
holographic models. The IR regime is dominated by the quadratic term and at T = the 
scalar field diverges logarithmically in the r = 0, near- horizon region. 

b) The potential behaves exponentially for small values of the radial coordinate r. Assuming 
that r = corresponds to — > oo, we have in this case 

V{cj,) ~ e^^ (2.12) 

where 6 is a positive constant. As we shall discuss later in Sect. llVCi this case produces a 
scale-covariant solution in the IR, corresponding to hyperscaling violation in the dual QFT. 
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c) The origin r = corresponds to an other extrema (a minimum) at cf> = (pi of the potential 
V{(l)). In this case the theory flows to a second AdS4 vacuum in the infrared. 

The IR regime of the EHT (|2.ip is also characterized by the IR behavior of the gauge coupling 
function Z. In particular, Z is crucial for determining the contribution of bulk degrees of freedom 
inside or outside the event horizon to the boundary charge density. This distinction is captured 
by the behavior of the electric flux in the IR 



$ = / Z{^)F , (2.13) 

where F is the dual Maxwell tensor. Using a terminology borrowed from condensed matter physics, 
the phase with $ = has been called cohesive and describes dual confined gauge invariant op- 
erators. The phase $ 7^ has been named fractionalized and describes a dual deconfined phase 
0jl23- -l-^^ this paper we will consider two choices for the gauge coupling function Z{(j)): (1) a 
minimal coupling, Z{(j>) = 1; (2) a coupling which behaves exponentially in the IR, Z ^ e"''^. 

Since in the following we shall make often use of the thermodynamical properties of the BB 
solutions, we find it convenient to summarize them here. The temperature T, entropy S and free 
energy F of the solutions (|2.2p are given by 



^^AJttO^ 5 = 47rVff2(r,0, F = M-TS, (2.14) 

47r 

where M is the total mass of the solution, V is the volume of the 2D sections of the spacetime and 
Th is the location of the outer event horizon. 



III. SPECTRUM OF EINSTEIN-SCALAR-MAXWELL ADS GRAVITY IN THE 

INFRARED REGION 

In this section we investigate general features of the mass spectrum of ESM-AdS gravity in the 
IR region. Assuming the existence of scalar-dressed BBs with AdS asymptotic behavior, the two 
basic questions in this context are about the existence and features of the T = extremal state 
and of the states near-extremality. We will treat separately the EM charged and uncharged cases. 
We will first consider the theory with zero charge density (Z ~ Y = in the action (|2.ip ). i.e. a 
vanishing Maxwell field (Einstein-scalar AdS gravity). Later, we will extend our considerations to 
the case of finite charge density. 



A. Einstein-Scalar AdS gravity 

A nontrivial point is the determination of the total mass M (i.e. the energy) of the BB solution. 
As discussed in Ref. J4j|, the usual definition of energy in AdS diverges whenever Oi ^ (with 
a divergent term proportional to r). This is because the backreaction of the scalar field causes 
certain metric components to fall off slower than usual. However, this divergent term is exactly 
canceled out if one considers that for Oi 7^ there is an additional scalar contribution to the 
surface terms which determines the mass. 

Using the Euclidean action formalism in the case m^ = —2/L'^ the total mass turns out to be 

m 

M ^Mg + ^ [O1O2 + P(Oi)] , (3.1) 

where Mq is the gravitational contribution to the mass, we have chosen the following boundary 
conditions for the scalar: O2 ~ 02{0i), and -P(Oi) = J 02{0i)d0i. 
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In the following we will need an expression for the mass when mg is in the range of values 
considered in this paper, —9/4 < rnj.L^ < —2. Furthermore, working with the parametrization of 
the metric given by Eq. ()2.2p , it is useful to express the total mass M in terms of the coefficient of 
the 1/r term in the r = oo expansion of the metric functions. To derive such an expression we use 
the Euclidean action formalism of Martinez et al. |4J|. Using the parametrization of the metric 
(12. 2 p the gravitational and scalar part of the variation of the boundary terms are given respectively 



SIg = —[{HH'8\-\'H5H) + 2\H{5H')]\^, 
T '* 

51^ = ^H^\cly'5<p\Z. (3.2) 

From the definition of the free energy F = M — TS, taking into account that /^Ir^ — 0, S — Iclrh 
and from F = —IT, it follows 

M = TS-TI^-T{I^ +1^). (3.3) 



To calculate the mass M p.3p we need the subleading terms in the r = oo expansion of the 
metric (|2.8p . By means of a translation of the radial coordinate r, the asymptotic expansion of the 
solution can be put in the general form: 



A = 


;Upr^-;;+o(.^- 


r2 ^ 


Lj t 


- 


%^%^Oir--^~^), 



(3.4) 

where p, q, a, /3, s are constants. Inserting this expansion in the field equations one gets (at the 
first and second subleading order) the following relations between the constants: 

R on A ^ AiQj A1A2Q1O2 ,„ ,. 

/3 = a = 2(l-Ai), l^ = g= 4^2(1 _2Ai) ' ' = ^I^' ^^'^^ 



Substituting Eq. (|3.4p into (|3.2[) and using p = q, we obtain 

-5/^ - -^ {Smo + 6ds - 26p{P - 1^+') , (3.6) 

SI^ - -^ (AiOi,50ir'9+i + A1O1SO2 + A2O2SO1) . (3.7) 

Notice that both the gravitational and the scalar contribution to the mass contain a term which 
diverges as r^~^^. Using Eq. p.Sp one easily finds that the two divergent terms cancel out in 
SI°° = SIq' + SI?' . Finally, we obtain the total mass of the solution 

M = -TI- = ^ (m„ + l^L-^ I dO^Wm + ^^^^^020^^ , (3.8) 

where we have parametrized the boundary conditions for the scalar in terms of the function Oi = 
1^(02)- It is also useful to split the total mass into the gravitational and scalar contributions Mq 
and M^, arising separately from the two terms in Eq. (|3.8|) : 

Mg ^ ^ (mo - ^^0102], M^ = -^[AiOi02 + (A2-Ai)P(Oi)], (3.9) 



where P(Oi) is defined as in Eq. p.ip . One can easily check that the previous equations reproduce 
correctly Eq. ((37T|) in the case m^ = -2/L'^, i.e. Ai = 1, A2 = 2. 

Let us now investigate general features of the mass spectrum of ES-AdS gravity in the IR region. 
In particular, assuming the existence of scalar-dressed BBs with AdS asymptotic behavior, we wish 
to characterize the features of the T — Q extremal state and of the near-extremal states. 

In the uncharged case a general, albeit implicit, form of the solution for the metric function A 
in a generic ES-AdS gravity theory has been derived in Ref. |34l |: 

X = H'-C,H'I^^, (3.10) 

where Ci is an integration constant. The equation above implies that if an extremal T = hairy BB 
solution exists, this must have Ci = 0, i.e. A = H^. We can prove this statement by the following 
argument. Differentiating the equation above and using Eqs. (|2.14p . we find the following relation 
between the temperature and the entropy density S of the solution: 

A'(r,,) (2XHH'-Ci),. 
j.^±SLJiI = 1 ^ZlA. (3.11) 

An S 

Therefore we get Ci = [2XHH']^ — ST. An extremal solution satisfies T = and A(r/j) = 0. 
Assuming that H and H' are finite at the horizon (to avoid curvature singularities), the existence 
of an extremal solution imposes Ci — 0, i.e. 

X = H^. (3.12) 

Note that this argument applies both when the entropy of the extremal solution is vanishing or 
when it is finite. 

Obviously, an extremal uncharged solution with AdS asymptotics (besides the trivial AdS vac- 
uum) may not exist. For the moment, we assume such a solution exists and derive a general and 
very important result. We shall prove that if such a solution exists it must have zero energy, i.e. 
must be degenerate with the AdS vacuum. 

To prove this statement, let us first show that every scalar-dressed solution with AdS asymptotics, 
which is characterized by A = H^ , requires necessarily conformal boundary conditions (j2.9p for the 
scalar field. The field equations (|2.3p - (|2.6p with p = and the metric (|2.2p are invariant under the 
scale transformation r — >■ fir, X — > /i^A, t — > /i^^t. In the extremal case, the asymptotic expansion 
(|3.4p implies that full solution (A = H^ and (f>) is invariant under this scale transformation if 
Oi^2 scale as follows: Oi -^ /i'^^Oi, O2 — > /i'^^02, which in turn implies the conformal boundary 
condition ((2^ . 

We can now calculate the mass p.Sp of the extremal solution, which has A — H^, hence mo = 
-2s. We get, 



L4 



(A2 - Ai)P(Oi) + Ai(l - ^A2)020i 



(3.13) 



where P{Oi) is defined as in Eq. (|3.ip . From Eq. (|3.13p it follows almost immediately that for 
conformal boundary conditions (|2.9p the mass M vanishes. 

This is an important and extremely nontrivial result. It means that in ES-AdS gravity with no 
EM field, if an extremal scalar-dressed BB solution exists then the AdS4 vacuum of the theory must 
necessarily be degenerate. Physically, this degeneration is a consequence of the fact that the scalar 
condensate gives a negative contribution to the energy. Therefore we may have configurations 
in which the positive gravitational energy is exactly canceled by the negative energy of the scalar 
condensate. This cancellation is a consequence of the conformal symmetry of the extremal solution; 
it necessarily occurs because the extremality condition A = H^ forces the conformal boundary 
conditions (|2.9p . 

It is also important to notice that the argument leading to the degeneracy of the T = ground 
state holds true also when a condition much weaker than Eq. (|3.12p is satisfied: 

X = H^ +0{r-^). (3.14) 
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In fact the mass p.8p and the scaling arguments leading to the conformal boundary conditions for 
the scalar field depend only on terms up to 0{r^^) and are completely insensitive to higher order 
terms in 1/r. 

Let us now consider near-extremal solutions. We assume that the theory allows for scalar-dressed 
BBs at finite temperature with AdS asymptotics. In the next section, we shall prove the existence 
of finite temperature solutions, by constructing AdS-BBs, numerically, for three classes of ES-AdS 
gravity models. 

The BB spectrum near-extremality can be investigated by considering the T — > limit of the 
finite T solutions. However, one can show that this T -^ limit is singular. In order to prove the 
statement we expand the fields in the near-horizon region, 

OO OO C30 

A = ^a„(r-nO", H = Y,bn{r~rhr, - ^ c„(r - r,,)" . (3.15) 

n— 1 n— n— 

At first order we get for 6o 7^ 

62 = -^c^, biai = --^V{co), '^i"=i=("^) ' (3-16) 

whereas the temperature of the dressed solutions, from Eq. p.I4p . becomes 

ai _ boVjco) 

Because in the case under consideration {V has a maximum) the potential V is limited from 
above {V{4>) < — 6/L^), the T ^ limit can only reached by letting bo — > 0. But on the other 
hand from the third equation in p.l6p it follows immediately that oi = is a singular point of the 
perturbative expansion p.l5[) unless {dV/d(j>)co = (corresponding to the Schwarzschild-anti de 
Sitter (SAdS) BB). Thus the T — >■ limit is a singular point of the perturbation theory. It should 
be stressed that this result has been derived by first considering the near-horizon limit, then taking 
T — >■ 0. In section flV Dl we will see what happens if the two limits are taken in the reversed order. 

Note that the above results are strictly true only if one considers AdS solutions with negative 
mass squared for the scalar field. If the scalar potential has a local minimum at = 0, then our 
argument above does not apply. This is for instance the case with the class of models studied in 
Ref. [33 which, however, turn out not to have BB solutions with AdS asymptotics. 

The singularity of the T — > limit in the near-horizon perturbation theory, indicates that the 
ground state (|3.12p is isolated, i.e. it cannot be reached as the T ^ limit of finite— T scalar BB 
solutions. This conclusion can be also inferred by reasoning on the r = 00 boundary conditions 
for the scalar field. We have previously shown that the symmetries of the field equations together 
with Eq. p.I2p force conformal boundary conditions ()2.9p for the scalar field. On the other hand, 
one can easily show that in the case of T finite, the field equations together with the conditions for 
the existence on an event horizon imply boundary conditions of the form (|2.10p , hence in general 
nonconformal boundary conditions. In fact, the field equations (|2.3p - (|2.6p have the following 
symmetries: 



r -> 


kr, 


t -> kt, 


L -^ kL, 


H^ 


kH 


r -> 


kr, 


X^k^ 


A, t^k^^t, 


Ao 


-)■ kAo 


H^kH, 


{x,y) 


-^ k-^{x,y), 






X^ 


fcA, 


t^ k' 


H, ij2 ^ k 


-'H\ 


, L^ 



(3.18) 

These symmetries can be used to fix all but one parameter in the perturbative expansion p.lSp . 
The solutions become in this way a one-parameter family of solutions. The near-horizon expansion 
(j3.15p depends on a single free parameter, which can be chosen to be r/;. For each value of r^, we 
can extract the two functions Oi(r/j) and 02(^/1)1 which define implicitly the boundary condition 
Oi = W{02). 
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It follows that in general the finite— T solution require boundary conditions for the scalar, which 
are different from the conformal ones required for the ground state p.l2[) . Therefore, the solution 
(|3.12p is generically isolated, i.e. it cannot be reached as the T ^ limit of finite— T scalar BB 
solutions. 

It should be stressed that the above feature is a key general feature of the BB solutions of 
AdS Einstein-scalar gravity which holds true also for the numerical solutions discussed in the next 
sections. If one assumes an analytic expansion close to the horizon, an asymptotically AdS behavior 
at infinity and if one requires the existence of hairy BB solutions, then the boundary conditions at 
infinity cannot be arbitrarily imposed but are determined by the field equations. These boundary 
conditions will have the form (|2.10p . with the function W determined by the dynamics. In the 
dual QFT the function W characterizes the scalar condensate. Thus, the particular form of the 
condensate is determined by the gravitational dynamics. It is obvious that this is true only in the 
case of pure Einstein-scalar gravity. For instance it does not hold for electrically charged solutions'^. 
In this latter case the near-horizon solution has always more than one free parameter, that can be 
fixed by prescribing some boundary condition for the scalar field. 

We can also compare the temperature of the dressed solution of radius rh with the temperature Tq 
of the SAdS BB with the same radius. We can use Eqs. (|3.18|) to set rh = L, 61 = ci = L^^, &o = 1, 
so that the only free parameter is cq = (t'(j'h) and the temperature becomes — SttT — LV{cq). We 
have therefore 

r - To = 87rL(l/(0) - y(co) ). (3.19) 

In the case under consideration, V{4>) has a local maximum at (j> ~ Q, so that V{Q) > V{co). 
Therefore, we obtain T > Tq for any finite temperature solution. That is, there exists a critical 
temperature given by the temperature of the SAdS BB: To = '^l2 such that scalar-dressed solutions 
of the same radius r/j only exist when T > Tq. 



B. Einstein-Scalar-MELXwell AdS gravity 

Let us now consider the EM charged case, i.e. a finite charge density in the dual QFT. In general, 
one expects that the finite charge density will remove the degeneracy of the T = extremal state 
we have found in the uncharged case. This can be shown explicitly. Indeed, when p 7^ 0, the field 
equations imply 

2 
' ' 2Ai/'^ + 2\HH" = H^\" , (3.20) 



zm 



which is solved by Eq. p.lO|) only when the charge is vanishing. In particular, A = H^ is not a 
solution of the equation above when p ^ 0. Moreover in the charged case Eq. p.lOp becomes |34| 



X = H' 



dr .2f,Jj_fj!_\ 



'-^^ljr^+p'Jd\H^j zm)\ 



(3.21) 



By using the same procedure leading to Eq. p.l2p . we get that the extremal solution in the EM 
charged case is attained for 






Ci-P' ^77H ■ (3-22) 



^ It does not hold also for black hole solutions of ES-AdS gravity, i.e. for solutions which spherical horizons. This 
is because in this case the field equations are not anymore invariant under the full set of transformations I I3.18I I . 
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Because Ci ^ 0, not even the weaker condition p.l4p is satisfied in the charged case. This implies 
that Eq. (|3.13p does not hold if p 7^ 0. In general, the mass of the extremal scalar-dressed solution 
will be different from the mass of the extremal RN-AdS solution, so that the degeneration of the 
r = ground state in the EM charged case is removed. 

Notice that in the charged case the T = solution is not necessarily forced to have conformal 
boundary conditions for the scalar field. In fact, the argument used for the uncharged case is based 
both on the relation A = H^ and on the scale symmetries of the field equations. Both do not hold 
anymore at finite charge density. Nevertheless, in this case the presence of an additional field (the 
EM potential ^^) allows to choose arbitrary boundary conditions for the scalar. As discussed in 
the previous section, the boundary conditions are not anymore imposed by the dynamics of the 
system as in the uncharged case. It follows that in the charged case the T = ground state is not 
anymore isolated but can be reached continuously as the T ~ limit of finite temperature scalar 
dressed BB solutions. 

For what concerns the BB spectrum near extremality, the results we have found in the uncharged 
case still hold in the charged case. In fact the first and the third equation in (|3.16p are not modified 
by the nonvanishing charge, whereas the second becomes aibi ~ —{Z^^{co)p'^ + 26o^(co))/45o. 
We obtain the temperature: 



An IGnbibf, 






Also here, we can compare the temperature of a scalar-dressed BB with that of the RN BB with 
the same charge p and radius r^,. One easily finds that Eq. p.l9p still holds for the charged case 
and that T > T^^ for any finite temperature solution, where the critical temperature T^^ is 
given by 

.^ _ 12rt-L'Z-\co)p' 

^° - le^^p^ ■ ^^-^^^ 

Scalar dressed EM charged solutions of the same radius r^ and charge of the RN-AdS solution 
exist only for T > T^^ . 

An important issue when dealin g w ith finite EM charge density is the characterization of the 
phase as fractionalized or cohesive [ia,l2^- For the generic theory (|2.1[) with F 7^ 0, this charac- 
terization will depend on the IR behavior of both Z{(f>) and Y{(j)). However, one can easily show 
that in the case of unbroken C/(l) symmetry, Y = 0, only the fractionalized phase may exist. In 
fact using Eq. (|2.6p into Eq. (|2.13p one easily finds <& ~ p. 

To summarize, the following interesting picture emerges for the IR spectrum of scalar-dressed 
BB solutions of ESM-AdS gravity with rri^p < ml < rn\p + XjL? . If a scalar-dressed, neutral, 
extremal solution exists at p = 0, it must necessarily be degenerate with the AdS vacuum. This 
is due to a precise cancellation of the contributions to the total energy from the gravitational and 
scalar part and, in turn, it is due to the conformal symmetries of the boundary theory. Moreover, 
the T — )■ limit of finite-T BB solution is singular and the ground state is isolated from the 
continuous part of the spectrum. 

When an EM charge is switched on, the degeneracy of the ground state is removed and the 
ground state can be reached continuously as the T — >■ limit of finite— T solutions. Scalar dressed 
uncharged (EM charged) solutions of the same radius r/j (and charge) of S-AdS (RN-AdS) solution 
exist only for T > Tq {T > T^^). Cohesive phases may exist only when the U{1) symmetry is 
broken. In the U{1) symmetry-preserving case only the fractionalized phases are allowed. 

Our results are fairly general and only assume the existence of scalar dressed solutions, which 
has to be investigated numerically. In the next two sections we will show that the picture above is 
realized for three wide classes of models with quadratic, quartic and exponential potentials V{(j)) 
and for two classes of gauge couplings [Z = \ and Z ^ e""^). Numerical computations confirms 
the degeneracy of the ground state in the uncharged case and the peculiarity of the T — >■ limit 
of finite-temperature scalar-dressed BB solutions. We will discuss separately the EM neutral and 
charged solutions. 
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IV. NEUTRAL SOLUTIONS 

A. Quadratic potential 

In this section we will construct numerical solutions of ES-AdS gravity models with the quadratic 
potential (|2.1ip and we shall check the validity of the general results of Section IIIII The case of a 
quadratic potential is the simplest possible choice and it is therefore our first example. Moreover, 
this is the usual choice for models describing holographic superconductors. We will come back to 
this point later in Sect. IV CI 

1. Extremal solutions 

The near-horizon behavior of the extremal solution of the mod el (12. lip with an EM field covari- 
antly coupled to a charged scalar field has been derived in Ref. [4^. The near-horizon, extremal 
solution of a pure Einstein-scalar gravity model (both the EM and the charge of the scalar field are 
zero) can be obtained as a particular case of the solution given in Ref. [45|. In the gravitational 
gauge used in Ref. 05|: ^^ have 

dp 
ds^ = -g{f)e~^de + — - + P{dx^ + dy^), (4.1) 

9{r) 

and with our normalization for the kinetic term of the scalar field, the solution reads 

ds'= ..f\ ^. +Pi-dt' + dx'+dy'), <^ = 2^/2(-lnf)i/2, go^-^. (4.2) 

gor^(— inr) 6 

The near-horizon, extremal solution (|4.2p can be written in the gauge (|2.2p by a suitable 
reparametrization of the radial coordinate. We get 

A = i?2 = e-^^, (b = ~j2^X{r), r=J—+ dte'^^, (4.3) 

V 50 J 

where the last equation defines implicitly the function X{r). We note that also in these coordinates 
the horizon is located at r = 0. 

The global, extremal, solution interpolating between the near-horizon behavior (j4.3p and the 
asymptotic AdS behavior p.4p has to be found numerically. We have integrated the field equations 
numerically for several values of m^. In all cases we have found A = i/^, which implies the conformal 
boundary condition (|2.9p . Indeed the total mass M of the scalar-dressed solution is zero. In Fig.[T] 
we show the profiles of the metric functions and the scalar field for m^ = — 2/L^. 

2. Finite-temperature solutions 

Let us now consider BB solutions of Einstein-scalar theory at finite temperature. Again we have 
to construct global solutions, which interpolate between the asymptotic AdS expansion given by 
Eq. (|3.4p and a near-horizon expansion as in Eq. p.lSp . 

We have constructed these solutions numerically, starting from the near-horizon solution above 
and integrating outwards to infinity, where the asymptotic behavior of the solution is AdS4. In 
Fig. [5] we show an example of the metric and scalar profiles and of the function 02(0i) in the 
case m^ = —2/L^^ . In the large Oi limit, our data are well fitted by O2 ^ — 0.57OJ, which 



^ In Fig. [2]and in all the figures we show in this paper all the dimensional quantities {Oi^2, F^J^, c,T) are normalized 
with appropriate powers of the AdS length L 
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r/L 



FIG. 1: The metric functions A = H^ (left) and the scalar field (j) (right) as functions of r/L, in the 
extremal case, for a quadratic potential with m^L^ = —2. 



is consistent with the conformal boundary condition (|2.9p . However, for small values of Oi the 
behavior reads O2 ~ — 0.36Oi and the global behavior interpolates between these two asymptotic 
regimes. Therefore, the function 02{0i) does not generically satisfy the conformal boundary 
conditions (|2.9I) . This is a general statement that we have verified also for different choices of the 
parameters and different models. This fact confirms that extremal solution are isolated from finite 
temperature solutions. 

As expected, solutions dressed with scalar hairs only exist above a certain critical tempera- 
ture Tq = 3/(47r) and a critical mass Mq which correspond to the temperature and mass of the 
Schwarzschild-AdS BB after a rescaling that sets rh = L = 1. This is shown in Fig. [31 where we 
present the total mass M of the solutions as a function of the temperature T. The absence of 
dressed solutions (irrespectively of the boundary conditions 02(0i)) for T < Tq confirms numeri- 
cally the existence of the critical temperature Tq. 
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FIG. 2: Left panel: metric and scalar profiles as functions of the (nonrescaled) coordinate r for a quadratic 
potential with irigL = —2. Right panel: the function O2 = 02(0i) for the same model. 
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FIG. 3: Total mass M as a function of the temperature T for a quadratic potential with TTigL 



B. Quartic potentials 



In this section we will check numerically the results of Section HIH and the validity of the picture 
that has emerged from our results, in the case of a theory with a potential 1^(0) having the behavior 
described as type c) in Section |lll i.e. a theory with a IR fixed point. As an example of such a 
theory we take the quartic potential, 



4j,4 



V{<P) = k'^cj, 



m 



6 



(4.4) 



This potential has the typical mexican hat form with a maximum at 



and a minimum at 



»12 



= ± 



2A2 



with V{(t>i) = -6/t 



= with V{0) 
= -mV(16A4) 



-6/i2, V"{0) = -m^ 

6/L^, V"{if)i) ~ 2171^. The potential is invariant under the discrete transformation — > —0, so 
that we will just consider </> > 0. The theory allows for two AdS4 vacua: an UV AdS at = 
(corresponding to r = oo), with AdS length L and with squared mass of the scalar given by — to^, 
and an IR AdSi at = 0i (corresponding to r = 0) with AdS length I and with squared mass of 
the scalar given by 2m^. Again, we focus on —9/4 < —m^L^ < —2. 



Extremal solutions 



A scalar-dressed, extremal solution of the kind discussed in the previous section would represent 
a flow between an UV AdS4 and an IR AdS4. Let us first investigate numerically the existence 
of such a solution. If it exists we know from the results of the previous section that it must have 
zero mass, i.e. it must be degenerate with the (UV) AdS vacuum. In order to construct such 
solution numerically we need its perturbative expansion in the UV (near r = oo) and in the IR 
(near horizon, r = 0) . The UV expansion is given by Eq. p.4p . For what concerns the near- horizon 
r = expansion, the field equations (|2.3p - (|2.5p give instead 



7^.3 



A^^-^.^ + 0(.^),i^4-^.^ + C^(. 



. = <^i + |r + 0(r2), (4.5) 



where 7 is an arbitrary constant. Moreover, Eq. (|2.5[) constrains the possible values of the param- 
eter m in Eq. (|4.4p . rn^ ~ 2/ P. Introducing a dimensionless parametrization for A in Eq. (j4.4p . 



A ^ = fcP, one finds that the restriction on m^ implies 



< fc < 



P _ k 



15 



(4.6) 



We have integrated the field equations numerically, starting from r ~ outwards to infinity. 
When (f> >Q, regular solutions only exist for 7 < 0. These solutions interpolate between the r = 00 
AdS behavior p.4p and the near horizon solution (|4.5p . 

In Fig. m we show the profiles of the metric functions and of the scalar for fc = 1. and the function 
02(01) (obtained by varying the free parameter 7) for selected values of fc. Again we have found 
that A = H^, which implies the conformal boundary condition (|2.9p and that the total mass M of 
the scalar-dressed solution is vanishing. 
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FIG. 4; Left panel: profiles of the metric functions and the scalar field, in the extremal case, for a quartic 
potential with fc = 1 and 7 — —1. Right panel: the function O-ziOi) for three different values of fc. 



2. Finite-temperature solutions 



Using the same method described in Sect. IIV A 21 we have constructed, numerically, dressed BB 
solutions at finite temperature for models with the potential (|4.4p . We have generated global BB 
solutions for m^L^ = —2 and for several values of the parameter A. These solutions interpolate 
between the near-horizon expansion (|3.15|) and the asymptotic AdS4 form. In Fig. [S] we show an 
example of the metric and scalar profiles and the function 02(Oi) for the case mlL^ = —2 and for 
some selected values of A. As it is clear from Fig.[5l the function 02(Oi) displays a universal linear 
behavior at small Oi, which already confirms that the boundary conditions are not conformal for 
any value of A. In addition, for larger values of Oi the slope of 02(0i) depends on the quartic 
coupling. 

In Fig. [6] we show the total mass of the solution as a function of the temperature for fixed 
horizon radius r^ = 1 and L = 1. As expected the dressed solutions exist only for T > Tq, 
confirming numerically the existence of the critical temperature Tq. It should be noticed that we 
have generated the numerical finite-temperature solutions for values of the parameters rri^ and A, 
which are different from those used to generate the extremal solutions. The reason for this choice 
is a numerical instability of the solutions for positive values of m^. 
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FIG. 5: Left panel: metric and scalar profiles as functions of the (nonrescaled) coordinate r for a quartic 
potential with m^L^ = —2 and A"* — 0.3. Right panel: the functions O2 ~ 02{0i) for different values of 
A. 
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FIG. 6: Total mass M as a function of the temperature T for a quartic potential and for selected values 
of the parameter A. 



C. Exponential potentials 

In this section we will investigate the case of a theory with a potential 1^(0) having the behavior 
described as type 6) in Section HIl i-e. the potential behaves exponentially ^ e''* for —> 00 
(corresponding to r = 0). 

We search for scalar-dressed BB solutions that smoothly interpolate between an asymptotic AdS 
spacetime and a near-horizon scale-covariant metric. In the dual QFT they correspond to a flow 
between an UV fixed point and hyperscaling violation in the IR. In general these interpolating 
solutions cannot be found analytically but have to be computed numerically. To be more concrete 
in the following we will focus on a class of models defined by the potential 



F(0) 



h-^V 



[cosh(60) + 36^ 



11. 



(4.7) 
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This potential is such that the mass of the scalar is independent from the parameter b, ml = —IjL?. 
Moreover, it contains as particular cases h = l/v3, 6=1 models emerging from string theory 
compactifications, for which analytical solutions are known |34l. |36|. 



1. Extremal solutions 

The leading near-horizon behavior of the extremal solutions can be captured by approximating 
the potential in the — > oo region with the exponential form V{(l)) = — (l/6^)e'"^. In this case the 
field equations ([231) -(ES]) give [H 

A = «o (^] , H = (^) , <^ = 00 - bwln (^^ , (4.8) 

e'"P»r1 2 

Notice that ao > requires w > 1/2. This restricts the parameter range to 1/2 < w < 2 
(0 < 6^ < 3). This ansatz provides an exact solution to the equations of motion with an exponential 
potential —{l/b'^)e^'^ but only the leading near-horizon, extremal, behavior of the solutions with 
V{(j)) generic. Solution (|4.8p is scale covariant, i.e. the metric transforms under rescaling with a 
definite weight: 

r ^ kr, {t, X, y) -^ fci-™(i, x, y), ds^ -^ k^'^'ds^. (4.9) 

The extremal solution (|4.8p contains an IR length-scale r_. However, in the case of neutral 
BB the scaling transformations (|4.9|) may change this scale. The metric part of the solution is 
scale-covariant whereas the leading logr term of the scalar is left invariant. The only parameter 
that flows when IR length-scale r_ is changed, is the constant mode (j)^) of the scalar. 

To reduce the number of independent parameters, we can exploit the symmetries of the field 
equations previously discussed [cf. Eqs. p.l8p ] to fix L = 1 and ^o = in Eq. (j4.8p . So we can 
start from the more simple ansatz containing only one free parameter r_ . 

Starting from this scaling behavior near the horizon and imposing an AdS behavior (|2.8p for the 
metric and the scalar field at infinity, we have integrated numerically the field equations with a 
potential given by Eq. (|4.7p , with different values of the parameter < 6 < \/i. We have found 
BB solutions with scalar hair, that interpolate between the near-horizon (|4.8p and asymptotic (|2.8p 
behavior. 

In Fig. [7] we show the metric functions and the scalar field of these extremal BBs for h = 1/2 
and the function P{Oi) (obtained by varying the free parameter r_) for different values of the 
parameter b. Also in this case we have checked numerically that A ~ H^ and that the conformal 
boundary conditions -P(Oi) ~ 0\ are satisfied. We have also explicitly checked that the mass of 
the extremal solution vanishes. 

For the two cases b = 1/a/3 and b = 1 the extremal solutions are known analytically [3J|. They 
are respectively given by 

X rr2 (?- + r_)5 3 ^/3 , 

\ = H- = r^, 0= -log 

L'^ 2 \r + r^ 



X = H' = -^r, c^ = - log (^-^- j , (4.10) 

where r_ is a constant. From solutions (|4.10p one can easily derive the function P{0\) defining 
the asymptotic boundary conditions for the scalar field. We have P = (2/\/3)Of for 6 = l/\/3 and 
P = 0\ioTb^ 1. 

In order to compare these analytical solutions with those obtained numerically, we need to 
eliminate a linear term in the asymptotic behavior of A(7'). Taking into account this translation, we 
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have checked expUcitly that our numerical solutions with b = l/VS and 6=1 and the numerical 
calculated functions P exactly reproduce the analytical results. In general, the proportionality 
factor / depends on the value of b. We observe that for 6 < 1 / is negative, for & = 1 / = 0, while 
for b > I f becomes positive, as shown in Fig. [71 
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FIG. 7: Left panel: profiles of the metric functions and the scalar field, in the extremal case, for the 
potential H4.7|) with b — 1/2 and r- — 1. Right panel: the function 02(0i) for three different values of b. 



2. Solutions at finite temperature 



Following the same method as the one used in the previous subsections, one can generate generic 
hairy BB solutions with AdS asymptotics at finite temperature, i.e. solutions interpolating between 
the near-horizon p.lSp and the AdS (|2.8p behavior. We have generated numerically these BB 
solutions and found, as in the case of quartic potential discussed above, that for every value of the 
parameter 6 in the allowed range, they exist only for r^ > 1. This implies the existence of a critical 
temperature Tq below which only the SAdS BB exists. 

A summary of our results is presented in Fig. |S1 which is qualitatively similar to Fig. [5] for the 
case of quartic scalar potential. Again we have verified that the function 02{0i) does not define 
conformal boundary conditions (|2.9p for the scalar, i.e. the extremal solutions are isolated from 
those at finite temperature. 



D. Perturbative solutions near-extremality 



In Sect, nil Al we have seen that the T — > limit of finite-temperature BB solution is singular 
and that the ground state p.l2p is isolated from the continuous part of the spectrum. A way to 
gain information about the behavior near-extremality is to consider separately the near-horizon 
and near-extremal expansion. In general the two limits do not commute. In this section we will 
perform this perturbative analysis for the potential (|4.7p . Similar results can be obtained for other 
classes of potentials. 

We look for perturbative solutions of the field equations (|2.3p - (|2.5p in the near-extremal, near- 
horizon regime. The near-extremal regime is obtained by expanding the metric functions A and H 
and (/) in power series of an extremality parameter m, with tti — ^ when the temperature T — )• 
(or the BB radius r^ — > 0). On the other hand the near- horizon regime is obtained by expanding 
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FIG. 8: Left panel: metric and scalar profiles as functions of the (nonrescaled) coordinate r for the 
potential (|4.7|) with b — 1/2. Right panel: the functions O2 = 02(0i) for different values of b. 



the metric functions and the scalar field in power series of r — r/j. Because in general the two limits 
TO — )■ and r — >■ r/i do not commute we have to consider separately the two cases. 



1. r -5> rh, Th ^0 
We first expand A and H and (f> in powers of m: 

00 00 CXD 

A(r) = ^A„(r)m", H (r) ^ Y. ^r^^^)'^'" ' (/.(r) = ^ (/.„(r)r 



(4.11) 



n=0 



For small BB radius rh << L (or equivalently small T, i.e. T << l/L) we can truncate in the 
perturbative expansion ()4.1ip to first order in m. At leading order we find that Aq, Hq and (J)q 
must satisfy the same field equations (|2.3p - (|2.5p . At subleading order we find instead 



H'^ = --{2Hocf>'ocl>i + {(b'o)^H,) 



{2XoHi + Hl\i)" = 4 [Ao (Hoffi)' + ^iH^H'^] 
(2Aoi/i + Hl\,)" = -2 fc/.ii/^rf^ ^ 2HoH,Vi<t>o] 



dViM , .,2^ d^viM 



o'Pi- 



(4.12) 
(4.13) 
(4.14) 

(4.15) 



A solution of Eqs. (|4.12p - (|4.15p can be obtained by setting (f>i = Hi = 0, so that they reduce to 
{HlM)" ^Q, {H^HoXi)'^0, (H^^'oXi)' ^ 0. (4.16) 

Equations p.3p - (|2.5p for the near-extremal leading order functions Ao,-ffo,0o can be now solved 
as a near-horizon expansion in powers of r. The leading term in this expansion being obviously 
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given by Eq. ((i?5|) . 

<^oM = -fet(^ln— + ^7n ( — ) ■ (4-17) 




n— ^ 

(n) 

For each order in the r-expansion we can then determine the corresponding term Aj^ for Ai by 
solving Eqs. (|4.16p . One could worry about compatibility of the three equations ()4.16p . However, 
one can easily realize that for H^ — cit\ the system (|4.16p is always solved by Ai = C2r~'+^ with 
ci_2 constant. This follows from the first equation in p.3|) . which implies Hq oc 1/r. The leading 
order in the near- horizon expansion involves w, ao, /3o)7o- The symmetry of the field equations 
(|2.3p - (|2.5p under a rescaling of H allows to fix /3o = 1, whereas as expected w,ao turn out to be 
given as in Eq. (|4.8p . At this order Eqs. (|4.f6p give in turn 

Ai°' (xr-'"+i. (4.18) 

At the n— th order in the near- horizon expansion we find X]^ ex r~"'~"+^. The form of the 
near-extremal solution is therefore given by, 

+ 0(m2), H = Ha, = 00, (4.19) 

where Aq, -ffoi 0o a-re given by Eqs. (j4.17p . Assuming m < in the previous equation, we find that 
at leading order the relation between m and r/j is m oc J"^""" . Notice that this is an expansion in 
1/?'. This means that terms with higher n give smaller contributions for r — > oo. 



2. rh ^ ,r ^ rn 

This limit has been already discussed in section Fill AI The expansion in powers of (r — rh) is 
given by Eq. (|3.15p and at leading order the field equations (|2.3p - (|2.5p give the relations (|3.16p 
involving the parameters ai.2, &o,i,2, co,i. At the next to leading order we have three more pa- 
rameters as, 63,02 and three more relations. We are therefore left with 4 independent parameters 
60, Co, ai, r/j. As previously discussed, the field equations have the symmetries (|3.18p so that r^ is 
the only independent parameter. In principle, one can now expand a„(r/j), 6„(r^), Cn{rh) in powers 
of rti, substitute in Eq. p.lSp and reorganize it as the power expansion in m given by Eq. (|4.1ip . 
By retaining only the linear terms in m one could then compare the result with Eq. (|4.19p . Unfor- 
tunately, this is a very cumbersome task. Indeed, terms of order 0{rh) are generated at any order 
in the near- horizon expansion p.lSp . The problem has to be solved numerically. Numerically, one 
can look for global solutions interpolating between the near-horizon behavior p.lSp with a given 
Th and the AdS asymptotic solution (|2.8p . There is no guarantee that the solutions obtained in 
this way match Eq. (|4.19p . This is because the two limits r — > r/j and r^ — > do not commute. 



3. Near-extremal numerical solutions 

We have generated numerically, for the case of the potential (|4.7p . the solutions interpolating 
between the AdS asymptotic behavior (|2.8p and the near-extremal regime given by Eq. (|4.19p . In 
Fig. [H] (left panel) we can see the profiles of the metric functions and the scalar field for b = l/v^ 
and rh = 10^^ (corresponding to m = —10^*). 
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FIG. 9: Left panel: interpolating solutions between AdS at infinity and the near-extreme regime given by 
Eq. (|4.19|l with h = l/\/3 and r^ ~ 10~^. Right panel: difference between the free energy density of the 
near-extremal BB solution and the free energy density of the SAdS BB, for two values of b. 

We see that although the T ^ limit of the near-horizon perturbation theory is singular and 
isolated, global solutions obtained interpolating the near-horizon, near-extremal behavior (|4.19p 
with AdS4 exist also for T > 0. This is a manifestation of the noncommutativity of the near-horizon 
and near extremal limit. From the point of view of perturbation theory, the T ~ singularity 
means that the perturbative series in m (|4.1ip do not converge and that solutions (|4.19p are only 
perturbative solutions valid for Vh << L. 

The hairy near-extremal solutions shown in Fig. IH] describe small thermal perturbations of the 
extremal solution, but they do not describe the small-T limit of finite-temperature solutions. These 
results confirm that the ground state solutions p.l2p are not smoothly connected to the finite-T 
solutions, because of the existence of the discontinuity. 

Perturbative solutions in the small scalar-field limit can be also constructed. These kind of 
solutions are described in the Appendix. 



E. Hyperscaling violation and critical exponents 



The extremal T = hairy solutions found in Sect. II V CI for the case of the potential (|4.7p describe 
a flow between the near-horizon regime (IR) hyperscaling violating regime and an asymptotic AdS 
fixed point. From a QFT point of view, this translates into a hyperscaling-violation phase in the 
IR and a scaling-preserving phase in the UV. We can characterize the holographic features of this 
flow by giving the scaling exponents in the conformal (AdS) phase and nonconformal (hyperscaling 
violating) phase. The IR behavior is dictated by Eqs. (|4.8p . The UV metric is instead that of 
AdS4. 

To describe hyperscaling violation in four dimensions it is useful to consider the following 
parametrization of the scale covariant metric: 



ds" 



-'-{-r-^^'-^'>dt^ + dx'- + dr^), 



(4.20) 



where 6 is the hyperscaling violation parameter and z the dynamic critical exponent. Under the 
followinE rescaline of the coordinates: 



t^ XH, 



• Ax,; 



22 

the metric (|4.20p transforms as: 

ds -^ X^/^ds. 
Moreover, this scaUng transformation determines the following scaling behavior for the free energy: 

Fr^T^^^. (4.21) 

By a simple redefinition of the radial coordinate and a rescaling of the coordinates, we can write 
the metric (|48|l in the form ([QO]) . We obtain: 

ds^ ^r^^{-dt^ + dxf + dr^). (4.22) 

Comparing Eq. (|4.22p with Eq. (|4.20[) . we can easily extract the parameters 9 and z of our solution: 

z = l, e^'^—^. (4.23) 

\ — w 

Notice that we are now using dimensionless coordinates r, t, Xi, so that the IR length-scale r_ drops 
out from the solution. While the value z = 1 of the dynamic critical exponent is largely expected 
for uncharged solutions, we see that 6'<Oforl<i(;<2 and 6* > 2 for 1/2 < w < 1, while 9 
diverges for w — I (recall that in our case 1/2 < w < 2). This is in agreement with the null energy 
conditions for the stress-energy tensor, which require, for z = 1 and in the general case oi d + 2 
dimensions, either 9 < or 9 > d. Trivially, the parameters of the UV AdS conformally invariant 
solution are z = 1, = 0. 

From Eq. (|4.2ip . substituting the (|4.23p . we get that the free energy scales as: 

F-TT^. (4.24) 

We see from Eq. (|4.24p that the exponent of T is negative for < w < 1 or. equivalently, when 
9 > 2. So in this case the free energy diverges for T — > and the corresponding phase is always 
unstable. 

F. Thermodynamics of the near-extremal solutions 

The hairy near-extremal solutions discussed above can be interpreted as small thermal fluctu- 
ations of extremal T = hairy BBs. The thermodynamical features of these BB solutions - in 
particular the free energy and the specific heat - will provide important information about the 
stability of the ground state. Properties such as the scaling exponents are determined by the be- 
havior of the system at the quantum critical point, namely by the T = scale-covariant extremal 
near- horizon solution (|4.8p . On the other hand the stability properties are global features and they 
must be investigated using the global T ^ solutions. 

By Eqs. (|2.14p . the temperature and the entropy density of the near-horizon, near-extremal 
solution (|4.19p are given at leading order by: 

T=^-^aorr\ S^ ^'"^^„ r^. (4.25) 

Notice that in these subsections we are using dimensionless coordinates, so that the IR length- 
scale r_ drops out from our formulae, as in Eq. ()4.22p . and we set L — \. Temperature and 
entropy density are therefore also dimensionless. 

The scaling exponent of the entropy becomes negative when 1/2 < w < 1 (corresponding to 
1 < 6^ < 3), implying a negative specific heat and the corresponding solutions are therefore 
unstable. This is in agreement with the results of the previous subsection concerning the scaling 
of the free energy for w < \. Moreover, in this case small values of the temperature correspond 
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to high values of the horizon r^ and of the parameter m, so that we cannot obtain near-extremal 
solutions (in the sense of small temperature solutions) with rh << 1, which is the range of validity 
of the perturbative solutions (|4.19|) . 

For what concerns the entropy density and free energy density /"o of the SAdS BB we have 

So = 4^^^ ^0 - - (y) T3. 

We have derived numerically the free energy of the numerical near-extremal solutions as a func- 
tion of the temperature, for T << 1 . In Fig. [S] (right panel) we show the behavior of the free 
energy density of the hairy BB solution compared with the free energy density of the SAdS BB for 
two selected values of the parameter b (both such that 1 < w < 2), and for small values T « 1 of 
the temperature. We observe that the scalar-dressed solutions are energetically disfavored against 
the SAdS BB. This result can also be verified analytically by comparing Fq with the free en- 
ergy density !F of the hairy near-extremal solution, which can be expressed as a function of the 
temperature using Eq. p.l4p . 

G. Thermodynamics of the finite temperature solutions 

We have also computed the free-energy F and the specific heat c of the finite-temperature 
numerical scalar-dressed solutions derived in the previous sections for the case of the quartic (|4.4p 
and exponential (|4.7p potential. The results are shown in Fig. (TU] where we plot Ai^/Fo and c as a 
function of the temperature, with Ai^ = F — Fq. The free energy F is always larger than that of 
the corresponding Schwarzschild-AdS BB at same temperature and the specific heat is negative. 
Hence, these solutions are energetically disfavored against the undressed ones. 



V. CHARGED SOLUTIONS 

In this section, we will extend the numerical results previously obtained for neutral BBs in ES- 
AdS gravity to the case of finite charge density, i.e to the case in which an EM field is present 
in the bulk. We will focus our attention on models with exponential (|4.7p or quadratic (|2.1ip 
potential. We will discuss separately the cases of: i) minimal gauge coupling Z ~ 1; ii) exponen- 
tial gauge coupling in the L'^(l)— symmetry preserving phase; iii) Minimal gauge coupling in the 
C/(l)— symmetry breaking phase {Z = I, Y ^ 0). 



A. Minimal gauge coupling 

In this section we will construct numerical BB solutions for the model (|2.ip with Z ^ 1,Y = 
and the potential (j4.7p . As usual we discuss separately extremal and finite temperature solutions. 



1. Extremal solutions 

Following the same approach as the one used for the case of electrically neutral solutions, we look 
for numerical scalar-dressed BB's solutions interpolating between an asymptotic AdS spacetime 
and a near-horizon scale-covariant metric. Also in this case, the near-horizon behavior can be 
captured by approximating the potential (|4.7p in the (/) ^ oo region with the exponential form 
V{4>) = — (l/&^)e'"^. The field equations (|2.3p - (|2.5p give the scale covariant solution, which in 
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FIG. 10: Free energy (left panels) and specific heat (right panels) of the dressed BBs as a function of the 
temperature. Top and bottom panels refer to the theory with a quartic potential with m^ — —2 and with 
the exponential potential (|4.7p . respectively. 



the dual QFT corresponds to hyperscaling violation: 



"°l^ ■ "H^ 



b ( T 

/)o-^(u' + 2)lnf — 



(5.1) 



w 



ao 



4/i 



262 



62 



62w(it; + 2)' 



2e'"^°(3u;-2) 
62(w + 2) ' 



where p is the charge density of the solution. The solution above, together with the condition 
ao > 0, restricts the parameter range to 2/3 < w < 2 (corresponding to < ^2 < 2). We can 
exploit the symmetries of the field equations to fix i = 1 and 0o = 0, leaving r_ the only free 
parameter. We immediately note an important feature of this solution: in the limit p — >■ it does 
not reduce to the near-horizon solution (|4.8p obtained in the electrically neutral case. This means 
that the uncharged solution (|4.8p and the electrically charged solutions (|5.ip represent two disjoint 
classes of solutions. 

As usual, starting from this near-horizon scaling and imposing an AdS behavior (j2.8p at infinity, 
we have integrated numerically the field equations for different values of the parameter 6, finding 
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numerical solutions only for b > 1/2. In Fig. [Tl]we show the fields for 6 = 1. As expected, here we 
find in general A ^ H^ , hence the mass of the solution is nonvanishing and the degeneracy with 
the AdS vacuum is removed. 




1 





1 

1 




1 

\ 
1 1 


- 0(r) 


\ \ 


- ■ Ao'(r) ; 






\\ 


- 




'^^^^^^^^^^^ 


~ -----___ ," 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 
r/L 



FIG. 11: Metric functions (left panel) and scalar and Maxwell field (right panel), in the extremal case, for 
the potential (|4.7[) with 6=1 and r_ = 1. 



2. Finite-temperature solutions 

At variance with the extremal case, the charge of finite temperature solutions is a free parameter 
and the uncharged case is obtained setting p = 0. Using a straightforward extension of the 
numerical integration previously discussed, we can construct finite-temperature solutions at fixed 
charge density p. Some examples are shown in Fig. [T^] for the potential (|4.7p with 6=1. In 
the left panel we show the radial profiles of the fields, in the central panel we show the function 
O2 = 02(01) for different values of p, and in the right panel we shown the difference between the 
free energy of the dressed solution and that of a RN BB with same radius and same charge. Notice 
that, as already stressed, in the charged case the boundary conditions can be arbitrarily chosen. 
In particular, one can also choose conformal boundary conditions of the form Oi = 0. However, in 
the case at hand we have found that such conditions do not allow for scalar-dressed BBs. 

Similarly to the uncharged case, these dressed solutions are always energetically disfavored with 
respect to the undressed ones. 



B. Nonminimal gauge coupling 

As an example of a model with nonminimal gauge coupling we consider here the model discussed 
in Ref. Q- The gauge coupling Z and the potential V are Q 



2'(0) = 2^0 cosh a(/ 



vw 



-21^0 cosh? 



(5.2) 



In the IR {(f> — > 00) both the gauge coupling Z and the potential V behave exponentially. The 
model therefore belongs to the wide class of EHTs that fiow to a hyperscaling violating phase in 
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FIG. 12: Left panel: profiles for the metric coefficients, scalar field and EM potential as functions of the 
(nonrescaled) coordinate r for the potential (|4.7p with b — 1 and different values of p. Central panel: the 
functions O2 ~ 02(0i) for different values of p and b = 1. Right panel: difference between the free energy 
of the dressed solution and that of a RN BB with same radius and same charge. 



the IR [6|-|8|, Il9l . I28fl . The extremal solution of the field equations in the near- horizon approximation 
is the scale-covariant metric Q: 



A 



ao 



4(a + b) 



H ^ 



w = 2 - Ach, 
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(5.3) 



4+ (a + 6)2' 

2Voe'"f''>ri p^ 

{w + 2h)iw + 2h-l)' Yq 



a + b' 



2VaS"^'>{2-2h-bS) 
(w + 2h) ■ 



In Ref. (7| r = global solutions interpolating between the near-horizon hyperscaling violating 
metric (|5.3p and the asymptotic AdS4 geometry have been constructed numerically. Furthermore, 
numerical finite-temperature solutions have been found and their properties have been discussed in 
detail. In particular, it has been shown that below a critical temperature Tc the system undergoes 
a phase transition. The scalar dressed BB solution becomes energetically preferred with respect 
to the RN BB. 

The results of Ref. Q fully confirm the general results of Sect. IIIII The finite charge density re- 
moves the degeneracy of the T = solution we have in the uncharged case. Comparing the charged 
solution (|5.3p with the neutral solution (|4.8p one easily realizes that although the IR behavior of 
the two solutions belongs to the same class (hyperscaling violating) the critical exponents change. 
Moreover, the nonminimal coupling between the scalar field and the Maxwell field is such that 
the energy of the extremal scalar-dressed solution is smaller than that of the RN-solution. This 
determines an IR quantum phase transition between the AdS2 x R^ near-horizon geometry of the 
RN BB and the near- horizon scale covariant geometry (|5.3p . In the dual QFT this corresponds to 
a phase transition between a conformal and a hyperscaling violating fractionalized phase. 

Because the thermodynamical properties of the system at small temperature are essentially 
determined by the T = quantum phase transition, this also explains why the hyperscaling 
violating phase is stable at small temperature, below Tc- 

The near-horizon solutions for the charged BB (|5.ip and (|5.3p depend on the same IR length- 
scale r_ as the neutral BB solution (|4.8p . However, in the charged case the scaling transformations 
under which the metric part of the solutions is scale-covariant, changes not only the constant mode 
of the scalar 0o but also the charge density p. Thus, changing the IR scale r_ corresponds to a 
flow of the charge density p. As noticed already in Ref. [2^ , this is an irrelevant deformation along 
the hyperscaling violating critical line. 
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It is also interesting to notice the different role played in the quantum phase transition by the 
finite charge density and the nonniininial gauge coupling. The finite charge density lifts the de- 
generacy of the T = vacuum and changes the values of the critical exponents of the hyperscaling 
violating solution, but it is by itself not enough to make the hyperscaling violating phase energet- 
ically competitive with respect to the conformal AdS2 x R^ phase. Indeed, in the case of minimal 
gauge couplings discussed in the previous subsection, the energy of extremal RN-BBs is lesser than 
the energy of the scalar-dressed T = solution. It is the non minimal coupling between the gauge 
and the scalar field that makes the extremal RN solution energetically disfavored with respect to 
the extremal scalar-dressed solution. 



1. Hyperscaling violation and critical exponents 

In the case of a potential behaving exponentially in the IR the near-horizon, extremal solutions 
are scale-covariant for both zero or finite charge density and for both minimal or nonminimal gauge 
couplings. On the other hand, the critical exponents are affected by switching on a finite charge 
density. In particular in the case of charged solutions we will always have z ^ I. 

In the minimal case, after a redefinition of the radial coordinate and a rescaling of the coordinates, 
the metric (|5.ip reads: 

ds^ = r^ ( -r^^^=^ dt'^ + dr'^ + dx] 



from which we can easily extract the critical parameters: 

2(2 -2w) 



= 4, z = 



w 



We note immediately that the hyperscaling violation exponent 6* is a (positive) constant, indepen- 
dent from the parameters of the potential. The range of w implies z < 1, which is in agreement 
with the NEC conditions. Indeed the latter impose, for this values of 9 and z, the conditions z > 2 
or z < 1. Moreover we note that for 1 < w < 2, z is negative. 

On the other hand, for 2/3 < u; < 1 (corresponding to < z < 1), the free energy scales with a 
negative exponent: 

which implies an instability of the corresponding phase and a negative specific heat. 

Finally, we consider the case of a nonminimal gauge coupling given by Eq. (|5.2p . The critical 
exponents can be read off from Eq. (j5.3p , after a reparametrization of the radial coordinate which 
puts the metric in the form (|4.22|) . We have 



4c 2c(2 - 2w) 



2c -I ' (2c-l)(2-w;)^ 



while the free energy scales as: 



(2c 



F - T^^ = T =(2—2) . (5.4) 



C. Symmetry-breaking phase 

As an example of a model having a C/(l)-symmetry-breaking phase we consider here the model 
discussed in Ref. [l|-|3, |45|, which gives the simplest realization of holographic superconductors. 
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The gauge coupling is minimal, while the potential V and the function F((/)) in the action (|2.ip 
are quadratic [45J: 

777 

z((/.) = i, 1/(0) = _^</,^ r(0) = gV, 

where g is the electric charge of the complex scalar field whose modulus is cf). 

The metric and scalar field associated to the T = solution of the field equations in the near- 
horizon approximation are given as in the neutral case discussed in Sect. IIV Al i.e. by Eq. (|4.2p . 
whereas the EM potential is Aq = (?!)or'^(- logf)^/^ with 2^ = -1 ± (1 - 48(j^/m^)^/2. Numeri- 
cal, extremal solutions interpolating between the near-horizon solution (|4.2p and AdS4 have been 
constructed for q^ > |777^|/6 in Ref. [45|. Numerical finite-temperature solutions have been also 
considered [l|-|3|. In particular, it is well known that below a critical temperature the supercon- 
ducting phase (corresponding in the bulk to the scalar-dressed BB solution) becomes energetically 
preferred. 

The results of Ref. [ll-y, |45| for the holographic superconductors fully confirm our general 
results of Sect. IIIII The finite charge density removes the degeneracy of the T = solution in 
the uncharged case. Moreover, the nonvanishing coupling function Y gives a mass to the U{1) 
gauge field and makes the extremal scalar-dressed solution energetically competitive with respect 
to the RN-solution. The system represents an IR quantum phase transition between the AdS2 x R^ 
near- horizon geometry of the RN BB and the near- horizon geometry (|4.2p . In the dual QFT this 
corresponds to the superconducting phase transition [iHsl, |45|, which occurs below the critical 
temperature. 

Similarly to the nonminimal case, also here the finite charge density and the nonvanishing 
function Y play a very different role. The finite charge density simply lifts the degeneracy of 
the T = vacuum we have in the uncharged case. But it is the coupling between the scalar field 
and the EM potential Aq that causes the superconducting phase transition to occur at the critical 
temperature. It is also interesting to notice that in this case the finite charge density does not 
change the metric (and scalar) part of the IR solution, which is determined by the near-horizon 
solution and it is described as in the EM neutral case by Eq. 



VI. CONCLUDING REMARKS 

In this paper we have discovered several interesting features of scalar condensates in EHTs, which 
may be relevant for understanding holographic quantum phase transitions. In particular, we have 
shown that for zero charge density the ground state for scalar-dressed, asymptotically AdS, BBs 
must be degenerate with the AdS vacuum, must be isolated from the finite-temperature branch of 
the spectrum and must satisfy conformal boundary conditions for the scalar field. This degeneracy 
is the consequence of a cancellation between a gravitational positive contribution to the energy 
and a negative contribution due to the scalar condensate. When the scalar BB is sourced by a 
pure scalar field with a potential behaving exponentially in the IR, a scale is generated in the IR. 

Switching on a finite charge density p for the scalar BB, the degeneracy of the ground state 
is removed, the ground state is not anymore isolated from the continuous part of the spectrum 
and the flow of the IR scale typical of hyperscaling violating geometries determines a flow of p. 
Depending on the gauge coupling between the bulk scalar and EM fields, the new ground state may 
be or may not be energetically preferred with respect to the extremal RN-AdS BB. We have also 
explicitly checked these features in the case of several charged and uncharged scalar BB solutions 
in theories with minimal, nonminimal and covariant gauge couplings. In the following subsections 
we will briefly discuss the consequences our results have for the dual QFT and for quantum phase 
transitions. 
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A. Dual QFT 

One striking feature of the uncharged scalar BB solutions we discussed is that the boundary 
conditions for the scalar field are either determined by the symmetries (for the ground state) or by 
the dynamics (for finite-temperature solutions). Because the only free function entering the model 
is the scalar potential ^(0), this means that the information about boundary conditions for the 
scalar field is entirely encoded in the symmetries of the field equations and in V. Since the scalar 
field drives the holographic renormalization group flow, this fact has some interesting consequence 
for the dual QFT. 

We have seen in Sect. Illll that in the case of zero charge density the ground state for the scalar 
BB must be characterized by conformal boundary conditions. Prom the point of view of the dual 
QFT this corresponds to a multi-trace deformation of the Lagrangian of the CFT. This is a relevant 
deformation, associated to a relevant operator, which will produce a renormalization-group flow 
from an UV CFT to an IR QFT. The nature of the IR QFT is entirely determined by the self- 
interaction potential, V{(j)). In the case of the quartic potential (|4.4p - which is characterized by two 
extrema - the IR QFT has the form of a further CFT. In the case of the exponential potential (|2.12p 
the IR QFT is characterized by hyperscaling violation. In the case of the quadratic potential (|2.1ip 
the characterization of the IR QFT is much less clear because of the absence of scaling symmetries. 

The characterization of the dual QFT at flnite temperature is much more involved. In this case 
we have generically nonconformal boundary conditions for the scalar fleld and the asymptotic AdS 
isometries are broken. Nonetheless, an asymptotic time-like killing vector exists and both the UV 
and the IR QFT should admit a description in terms of multi-trace deformations of a CFT. 

On the other hand we have shown that the ground state and flnite— T states are not continuously 
connected. This means that we are dealing here with two different disjoint sets of theories. 

This picture changes completely when one adds a flnite charge density. Now the boundary 
conditions for the scalar fleld can be arbitrarily chosen, for instance in the form of the usual 
conformal Neumann or Dirichlet boundary conditions. Thus, in the case of finite charge, we have 
the usual description borrowed from the AdS/CFT correspondence with single trace operators dual 
to the scalar field. 



B. Scalar condensate and quantum criticality 

The results of this paper improve our understanding of quantum critical points in EHTs. In 
particular they shed light on the phase structure of these critical points proposed in Ref. [23| and 
on their stability. 

The degeneracy of the ground state for uncharged BBs simply means that at zero charge density 
the hyperscaling violating critical point (or line) and the hyperscaling-preserving critical point have 
the same energy. The potential V for the scalar field determines completely the scaling symmetry 
and the critical exponents of the hyperscaling violating critical point. The renormalization group 
fiow from the UV conformal fixed point into the IR introduces an emergent IR scale. Changing 
this IR scale produces a flow of the constant mode of the scalar field. As already noted in Ref. |28j . 
the presence of this arbitrary scale implies that hyperscaling violating critical points appear as 
critical lines rather than critical points. On the other hand, for scalar-dressed BBs the ground 
state is isolated from the finite— T part of the spectrum and the states at finite temperature are 
always energetically disfavored with respect to the SAdS BB. Thus, at zero charge density there 
is no phase transition between the hyperscaling preserving phase and the hyperscaling violating 
phase. 

Considering charged scalar BBs, i.e. introducing a finite charge density p in the dual QFT, 
generates several effects. First of all the degeneracy of the ground state is lifted and the ground 
state is not anymore isolated from the T > continuous branch of the spectrum. The change of 
the IR scale typical of hyperscaling violating critical lines now also produces a flow of the charge 
density p. Although the critical exponents are modified by the presence of a finite charge density 
(for instance the dynamical critical exponent z becomes ^ 1) the scaling symmetries characterizing 
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the critical point are very similar to those we have in the case of p = 0. The similarity between the 
ground state geometries in the p ^ Q and p = case is even more striking in the case of covariant 
gauge coupling (the case of holographic superconductors). In this latter case the metric and scalar 
part of the near-horizon solution is exactly the same for p ^ and p = 0. 

The stability of the hyperscaling violating critical line is a far more involved question. It turns 
out that it depends crucially on the coupling between the scalar condensate and the EM field, 
i.e. on the two coupling functions 2'(0) and Y{(t)) in the action p.l|) . In all cases that we have 
considered with a minimal gauge coupling Z = 1, and in absence of U{1) symmetry breaking 
(y = 0), the hyperscaling preserving phase is always energetically preferred with respect to the 
hyperscaling violating one. In this case, an IR phase transition between hyperscaling-preserving 
phase and the hyperscaling violating phase does not occur. 

Conversely, in the two cases of a nonminimal gauge coupling behaving exponentially in the 
IR {Z ~ e""^, y = 0) and covariant gauge coupling {Z = 1, Y ^ 0^), the hyperscaling violat- 
ing phase is energetically preferred. This gives, respectively, the IR phase transitions between 
the hyperscaling preserving phase and the hyperscalingviolating phase found in Ref. [7| and the 
well-known superconducting phase transition of Ref. [l|-y, \^. On the other hand, considering 
charged BBs at finite temperature, the critical temperature of the phase transition between the 
hyperscaling-preserving/hyperscaling violating phases is settled by the charge density p |7|, i.e by 
the IR emergent scale typical of the hyperscaling violating critical line. 

Summarizing, our results strongly indicate that for EHTs described by (|2.ip the three coupling 
functions V{(j)), Z{(j)), Y{(j)) determine different features of holographic quantum critical points. 
The self-interaction potential V{(j)) determines the scaling symmetries but not the stability of 
hyperscaling violating phases. Conversely Z and Y are crucial in determining the stability, the 
breaking of the U{1) symmetry and the characterization as fractionalized or cohesive of the hyper- 
scaling violating phase. 
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Appendix A: Uncharged perturbative solutions in the small scalar-field limit 

In the neutral case, it is possible to construct analytical BB solutions in the small scalar-field 
limit perturbatively, i.e. expanding the solution as follows 

r^ Af 
Kr) = II -^ + ^'^2(0, (Al) 

H{r) = r + e^i/aW, (A2) 

</>(r) = 6</>i(r). (A3) 

where e is a book-keeping parameter of the expansion. The solution for the scalar field can be 
obtained by solving the scalar equation at first order. The regular solution, can then be inserted 
into the Einstein equations that, to second order, can be solved for A2 and H2. 

Let us start with the T ~ AdS4 vacuum, i.e. we set iVf = in the equations above. To second 
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MO = -. + [-Tr^-TTh2+^^ + -]^' + 0{^'). (A4) 



order in the scalar field, the solution reads 



L2 ■ V 4L2 6L2r2 + 7,2 r 



^(0 = r + (-^ - ^ - ^ + ^^ + -^2) e^ + 0(.^) , (A5) 

<^M = ^f?^ + 5')+0(^')^ (A6) 



where Ci are integration constants. This is a solution for the classes of potentials presented in the 
main text. Although not presented, the solutions can be obtained in closed form at least to fourth 
order. The constant C\ can be set to zero by performing a coordinate translation such that the 
asymptotic form of the metric reads as in Eq. p.4p with Gi = —mQJI being related to the metric 
contribution to the gravitational mass and after a rescaling H — > -ff/(l + £^^2), which can be 
performed by rescaling the the transverse coordinates. Interestingly, there exists an event horizon, 
so the solution represents a BB endowed with a scalar field. Let us consider two cases separately: 
O2 = and O2 = 02(01) (without loss of generality, we assume O2 > 0). For the latter case, the 
horizon is located at 

Vg|. ^ 3V^ (40fO| 3,„g) ^.,, 
61/4 4^0^ 23/432O2'/' 

and, to first order, the temperature of the solution is 



T=^^^. (A8) 

On the other hand, if O2 — 0, the horizon and the temperature read 

_ 3my'£2/3 04g2 

^ - 21/34^ + ^^^ ^^^^^ 

In general, these solutions describe a BB whose horizon shrinks to zero in the O^ — >■ limit. The 
total mass of the BB is given by Eq. p.ip and it coincides with mo when O2 = 0. It is interesting 
to compare the free energy of this solution with that of a SAdS BB at the same temperature. 
When O2 7^ 0, we obtain 

^-^O = ^^g^ + O(a, (AH) 

so that F > Fq for any O2 7^ and the dressed solution is always energetically disfavored. Note 
that this result is valid for any boundary condition O2 = 02(0i) 7^ and for any scalar potential 
whose expansion reads V ~ —6/L^ — c/P /L"^. On the other hand, if O2 = 0, F = Fq to second 
order in Oi , so that the two solutions are degenerate. 

Finally, we can adopt the same technique to construct perturbative solutions of the SAdS BB 
at finite temperature. At first order, the general solution of the scalar field equation reads 

^1 = «P„i/3 [r'/iLHl) - 1] + /3Q_i/3 [r^/{LHl) - l] , (A12) 

where P„ and Qn are Legendre functions of order n and a and /3 are integration constants. Imposing 
regularity at the horizon rh = (2L2j\,/)i/3 requires /3 = 0. In principle, this solution can be inserted 
in the Einstein equation in order to obtain two equations for H2{r) and \2{t). Unfortunately, these 
equations do not appear to be solved in closed form. 
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